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ABSTRACT 


An  algorithm  is  described  for  reconfiguring  a  2-dimensional  VLSI  array  on  a  silicon 
wafer  that  has  some  faulty  cells.  The  functional  cells  of  the  array  are  intercon¬ 
nected  in  order  to  simulate  a  fault-free  array  of  smaller  size,  where  the  interconnec¬ 
tion  wires  are  routed  inside  horizontal  and  vertical  channels,  according  to  the 
Manhattan  model.  The  concept  of  simulation  distance  is  introduced,  and  it  is  shown 
to  be  related  to  the  length  of  the  longest  interconnection  wire.  The  algorithm  makes 
use  of  network  flow  techniques  in  order  to  find  a  wiring  with  minimum  simulation 
distance.  This  results  in  a  practical  heuristic  for  minimizing  the  maximum  wire 
length.  The  complexity  and  performance  of  this  algorithm  are  also  discussed  in  the 
paper. 


KEYWORDS:  network  flow.  VLSI  array,  fault  tolerance. 
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1.  INTRODUCTION 


The  technique  of  wafer  scale  integration  for  VLSI  circuits  has  received  considerable  attention 
in  recent  years.  The  basic  idea  is  to  assemble  a  large  system  of  processors,  or  cells,  on  a  single  sili¬ 
con  wafer  so  that  the  chip  packaging  costs  are  cut  off.  Due  to  the  physical  and  technological  limits 
of  the  integration  process,  some  cells  of  the  wafer  can  be  defective,  or  "dead".  Therefore,  the  prob¬ 
lem  arises  of  reconfiguring  the  interconnection  network  using  the  "live"  cells. 

The  reconfiguration  of  one  and  two-dimensional  arrays  of  cells,  typical  of  VLSI  systolic 
architectures,  has  been  investigated  in  [4.  2],  In  both  papers,  a  probabilistic  model  of  cell  failure  is 
adopted  and  algorithms  for  minimizing  the  maximum  wire  length  in  the  reconstructed  array  are 
given.  Channel  width  and  area  penalties  are  also  considered. 

In  this  paper  we  present  a  new  approach  to  the  problem  of  reconfiguring  two-dimensional 
arrays,  based  on  network  flow  techniques.  Section  2  illustrates  the  basic  ideas  on  which  our 
approach  is  based.  In  Section  3.  we  give  an  algorithm  for  the  reconstruction  of  two-dimensional 
arrays.  Section  4  contains  some  conclusive  remarks'. 


2.  PRELIMINARIES 

We  assume  the  same  model  for  the  wafer  scale  system  as  in  [4].  The  cells  are  positioned  in  a 
v fn  X'fn  array.  The  strips  between  two  rows  or  columns  of  cells  are  called  channels.  Each  chan¬ 
nel  contains  a  fixed  number  of  tracks  on  which  the  interconnection  wires  are  routed.  Suppose  that 
m  cells  in  the  array  are  dead.  We  investigate  the  problem  of  interconnecting  the  remaining  live 
veils  into  a  square  array  of  size  | 'Jn  —  m  |x|\/rc  —  m  |.  with  the  goal  of  minimizing  the  maximum 
wire  length  in  the  reconstructed  array.  See  an  example  in  figs.  1  and  2.  taken  from  [4] 

Consider  the  distribution  of  faults  shown  in  fig.  3.  It  is  easy  to  see  that  the  wiring  ol  fig.  4 
provides  an  optimal  reconstruction.  This  observation  suggests  the  following  reconfiguration  stra- 
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Figure  2  Example  of  reconfiguration  oi  the  arrav  in  fig.  1.  so 
that  the  live  cells  form  a  square  array  ot  3 b  ceils. 


I  1  )  Select  an  intermediate  array  whose  dead  cell  distribution  consists  of  exactly  \fn  —  Jn  —  m 
rows  and  •Jn  —  Jn  —  m  columns. 


TTTT 


(2)  Interconnect  the  live  cells  of  the  given  array  in  order  to  construct  the  intermediate  array 


(3)  Route  the  intermediate  array  by  connecting  each  live  cell  to  the  closest  live  cells  on  the  same 
row  or  column. 
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Figure  3  An  array  of  64  cells  whose  dead  cells  are  arranged 
into  two  rows  and  two  columns. 


Figure  4  Optimal  reconfiguration  of  the  array  in  fig  3 
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The  second  step  can  be  viewed  as  a  simulation  of  the  intermediate  array  by  the  original  one. 
Namely,  each  live  cell  of  the  intermediate  array  is  simulated  by  a  distinct  live  cell  of  the  original 
array. 

We  introduce  now  the  basic  terminology  on  network  flows,  to  be  used  in  the  next  sections. 

A  (flow)  network  is  a  6-tuple  N  =  (V  .A  ,b  ,c  .s  .t  )  where: 

( 1 )  U  is  a  set  of  nodes'. 

(2)  A  QV  xU  is  a  set  of  directed  arcs  (the  digraph  with  vertex  set  V  and  arc  set  A  is  called  under¬ 
lying  digraph  of  AO; 

(3 )  b  1’  [J  A  — 1 >  N  associates  to  all  nodes  and  arcs  a  nonnegative  integer  capacity', 

(4)  c  .A  -*  N  associates  to  each  arc  a  nonnegative  integer  cost'. 

(5)  s  and  t  are  two  designated  nodes  called  the  source  and  the  sink,  respectively. 

A  flow  for  A'  is  a  function  f  :A  -*R  that  satisfies  the  following  conditions: 

( 1 )  0  ^  j  (u  ,v  )  ^  b  (u  ,v  )  for  all  (u  ,v  )  €  A  ; 

(2 )  £  /  (u  .v  )  =  £  /  (v  ,w  )  ^  b  (v  )  for  all  v  €  V  —  { s  ,t  }. 

«■(«.'  1 1  A  v  •  (v  >•  ) €  A 

The  cost  of  the  flow  /  is  the  quantity: 

COST( /  )  =  £  c(u.v)  f  (u.v). 

(«  ,v  )  €  A 

The  value  of  the  flow  /  is  given  by: 

VAI.UEf/  )=  £  I  (s  x  ). 

I  '  N  1  '  t  \ 

A  minimum  cost  flow  of  value  <b  lor  N  is  a  flow  /  for  A'  such  that  VAl.L  lil  /  )  =  d>  and  COS  H  /  )  is 


minimum. 
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3.  THE  NETWORK  FLOW  TECHNIQUE 

We  represent  the  2-dimensional  array  with  an  undirected  graph  G  =  (V  £  ).  where: 

V  =  {(i  ,j  )  I  1  ^  i  .j  ^  \fn  ),  and 

E  ~  {((t  .j  )Xk  .1 ))  I  it  -k  I  +  \j  -l  I  =  1  }■ 

Let  dG  (u  ,v  )  denote  the  distance  between  vertices  u  and  v  in  the  graph  G.  The  dead  cells  are  in  a 
subset  of  V.  denoted  by  D. 

In  the  rest  of  this  section,  we  assume  that  the  intermediate  array  is  given  and  we  denote  by  T 
the  set  of  its  dead  cells,  which  we  call  target.  Note  that  I  T  I  ^  I  D  I .  The  simulation  between  the 
cells  of  the  original  and  intermediate  array  is  expressed  by  a  function  cr  mapping  vertices  into  ver¬ 
tices  where  cr(u  )  =  v  iff  cell  v  of  the  intermediate  array  is  simulated  by  cell  u  of  the  original  array. 

The  function  cr  must  satisfy  the  following  properties: 

( 1 )  no  cell  is  simulated  by  more  than  one  cell: 

(2)  only  the  live  cells  can  perform  simulations: 

(  3  >  no  dead  cell  in  the  target  is  simulated: 

( 4  )  all  cells  outside  the  target  have  to  be  simulated. 

More  formally,  this  can  be  expressed  by  defining  cr  as  a  one-to-one  function  (property  1) 
mapping  the  set  V  —  D  (property  2)  into  the  set  V'  —  (/)  f)r)  (property  3).  such  that,  for  all 
v  €  V  —  T  .  there  exists  a  vertex  u  such  that  cr(u  )  =  v  (property  4).  We  associate  to  cr  a  directed 
graph  (;  ,  =  1  \  ' n .A  „ ) .  where  \',T  =  V  —  (D  T  )  and  (u  ,v  )  €  A  ,T  w  henever  cr(u  )  =  v  . 

Lemma  1  A  function  <7  is  a  simulation  mapping  if  and  only  if  the  graph  G,r  consists  of  disjoint 
paths  that  are  either  cycles  of  (me  or  more  live  cells  or  chains  starting  f  rom  a  live  ceil  in  the  target, 
ending  at  a  dead  cell  outside  the  target,  and  having  the  intermediate  cells  all  lice.  Furthermore,  the 
number  of  these  chains  is  exactly  1  D  P|  (U  —  T  )  • 
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Proof:  Assume  that  O’  is  a  simulation  mapping.  Since  <7  is  injective,  each  vertex  in  V'^  has  at  most 
one  ingoing  and  one  outgoing  arc.  This  implies  that  G„  consists  of  disjoint  cycles  and/or  chains. 
Because  a  is  defined  on  the  set  of  live  cells,  then  each  cycle  in  G  „  contains  only  live  cells  and  every 
chain  in  Ga  consists  of  all  live  cells,  but  the  last  one.  which  is  a  dead  cell  outside  the  target. 
Furthermore,  all  chains  start  at  a  live  cell  inside  the  target,  since  each  vertex  in  V  —  T  has  an 
inverse  mapping.  Finally,  each  dead  cell  outside  the  target  has  an  ingoing  arc.  so  that  the  number 
of  chains  is  exactly  I  D  P)  (V  —T)  I .  This  completes  the  proof  of  the  Only-If  part.  The  proof  of 
the  If  part  readily  follows.  □ 

Lei  dr,  (u  ,v  )  denote  the  distance  between  vertices  u  and  v  in  the  graph  G.  i.e.  the  minimum 
number  of  edges  in  any  path  from  u  to  v.  The  maximum  wire  length  l„  in  the  interconnection  of  the 
intermediate  array  generated  by  cr  is  given  by: 

l„  =  max  jdG  (u  ,t>  )  i  (cr(u  ).cr(v  ))€£}. 

We  introduce  now  the  concept  of  simulation  distance  d defined  by: 

d „  -  max  {dG  (u  ,cr(u  ))  I  u  €  VCTj. 

which  is  related  to  the  maximum  wire  length  as  stated  in  the  following  proposition. 

Proposition  1  Let  cr  be  a  simulation  mapping,  then  1<T  ^  2d „  +  1. 

Proof:  A  consequence  of  the  inequality:  dG  (u  ,v  )  ^  dG  (u  ,w  )  +  dG  (w  ,v  ).  □ 

Given  a  positive  integer  k,  an  auxiliary  graph  can  be  used  to  check  w  hether  d  „  ^  k  .  Let 

G  ’  =U'  ..1 )  be  the  directed  graph  defined  as  follows. 

(1 )  (’  =  \  (J  Ij  ,(  |.  where  >  and  t  are  two  new  distinct  vertices. 

(2)  A  1  =,L  (J  A  f  '  U  A.  .  where: 

■  L  =  l(.»  M  )  u  €  T  n  (V  -D)\. 

A,:,1  =  Hu  ,v  )  €  xV  I  u  €  \  -  D  .  V  €  \  -  (D  f)  T  ).  and  d(!  (u  ,v  )  ^  k  \.  and 


>■1 

.J 
-  ■  1 
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A,  =  |(«  ;  )  I  u  €  D  f)('/  -D). 

The  length  of  arc  (u  .v  )  in  A,^  1  is  defined  as  the  distance  between  u  and  v  in  G.  Fig  5  shows 
the  digraph  G'1'  corresponding  to  the  array  of  fig.  1  and  to  the  target  of  fig.  3. 

Proposition  2  Let  G.  D  and  7  be  an  instance  of  the  simulation  problem.  Then  there  exists  a  simu¬ 
lation  function  cr  with  simulation  distance  d  c  =  k  if  and  only  if  the  digraph  Ga)  contains 
\  D  P)  ( V  —7)1  disjoint  directed  paths  from  s  to  t. 

Proof:  Follows  from  the  definition  of  G(t  >  and  Lemma  1.  □ 


In  order  to  minimize  the  global  wire  length  in  the  reconstruction  of  the  intermediate  array, 
we  look  for  a  set  of  disjoint  paths  in  Gu  1  with  minimum  global  length.  Moreover,  it  is  convenient 
to  choose  the  arcs  in  each  path  as  short  as  possible,  so  to  decrease  the  probability  that  la  attains  the 
bound  2d  „  +  1.  For  this  purpose,  we  can  associate  a  cost  c  (u  ,v  )  =  x(dG  (u  .v  ))  to  each  arc  (u  ,v  )  in 


Figure  5  The  digraph  G  1  corresponding  to  the  array  ol  tig 
1  and  to  the  target  of  fig  3.  The  dead  cells  outside  the  target 
are  represented  by  a  circle.  The  live  cells  inside  the  target 
are  represented  by  a  triangle.  Only  the  arcs  in  .L„‘  are 
shown,  where  an  arc  without  arrow  stands  tor  a  pair  ol 
svmmetric  arcs. 
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A  '•*  depending  on  the  length  of  the  arc,  so  that  the  following  conditions  are  satisfied: 

(1)  *w,+  ■  +dp  )  >  x(d  ,)  +  +  *(</.,  ).  lor  all  d,  €  N.  i  =  1 . p  . 

(2)  \(cj)+  •  +  x(cy  )  <  l)  +  +  *( dp  )  implies 

c  [4-  +  cq  <  d  j  +  •  +  dp  .  for  all  c,  .d j  €  N.  i  =  1 . q  .  y  =  1 . p  . 

Both  conditions  are  fulfilled  by  the  cost  function  x(d  )  =  di+(.  provided  e  is  a  sufficiently 
small  positive  constant. 

Now,  we  present  the  algorithm  SIM  that  finds  a  simulation  function  cr  with  minimum  simu¬ 
lation  distance,  given  a  positive  integer  n.  and  the  sets  D  and  T  of  dead  and  target  cells,  respec¬ 
tively. 

Algorithm  SIM 

Input:  positive  integer  n.  sets  D  and  T  of  dead  and  target  cells,  respectively: 

Output,  simulation  function  cr  with  minimum  simulation  distance. 

begin 

6  :=  I  D  f|  (v  ~T  )l : 

Bl’ILDNET  ( 1  .A’'1 '); 
if  MAX  FLOW  (A(1))  ^  0 
then  begin 

/  :=  M I XCOSTF LOW 
SIOMA  ( ,V  1  .  /  .cr  ! 

end 

else  begin 

a  :=  1 : 

b  :=  2 \ fn  —  2 : 
k  :=  | (a  +b  )/2|: 


BUILD  NET  (*  JVU  >): 
if  MAXFLOW  (Nu))  2  <f> 


then  b  \-k 
else  a  k  ; 

*  :=  |(a  +  6  )/2| 

until  a  =  b: 

/  :=  M  INCOST  FLOW  (Nik  ).<3b); 

S/GAM  (7V(*  >./  ,cr) 

end 

end 

procedure  BUILDNET  (k  .Nu  )); 
begin 

Build  the  auxiliary  flow  network  Nu  !  =  ((/  ..4  u  \b  ,c  ,.v  .f  ) 
as  follows: 

the  underlying  digraph  of  Na)  is  Gu  ): 

each  vertex  and  arc  has  unit  capacity: 

each  arc  ( u  .v  )  in  A„\k  ’  has  cost  c  ( u  .v  )  =  (dG  ( u  ,v  ))‘  + 

for  all  the  remaining  arcs  c  (u  ,v  )  =  0. 

end 

procedure  SIGMA  ( N“‘.j  .cr ): 

begin 

for  each  arc  {u  .v  )  €  do 
if  /  (u.v)=  1 


then  <t(u  )  :=  v ; 


function  MAX  FLOW  (N  ); 


begin 

Return  the  maximum  flow  value  for  the  networkA/ 

end: 

function  MINCOSTFLOW  {N  .</>): 
begin 

Return  a  minimum  cost  flow  of  value  <t>  for  the  networkA/. 

end: 

Figs.  6  and  7  show  the  application  of  algorithni  SIM  to  the  array  of  fig.  1. 


Figure  6  Minimum  cost  flow  in  the  network  .V  : 
corresponding  to  the  arrav  of  iig  1  and  to  the  target  ol  fig  3 
Only  the  arcs  in  A,,'.-  with  nonzero  (unit)  flow  are  shown 
\ote  that  k  -  2  is  the  minimum  value  tor  which  A/1*  ad 
mits  a  feasible  flow  of  value  i  D  P|  ( V  —  7  ) !  =  1 5. 


Figure  7  Reconfiguration  of  the  array  in  fig.  1  produced  by 
the  algorithm  SIM. 


Proposition  3  Algorithm  SIM  correctly  computes  a  simulation  function  <r  with  minimum  simula¬ 
tion  distance  in  O  (n3)  time. 

Proof:  The  correctness  of  the  algorithm  follows  from  the  fact  that  any  minimum  cost  flow  for 
,V  *  of  value  0=l£>f'|(V— r)l  consists  of  <f>  disjoint  paths  from  r  to  t  [5]. 

At  most  O(logn)  maximum  flow  computations  are  performed,  each  taking  time  Oin-3)  [l. 
cap.  6].  The  minimum  cost  flow  computation  is  executed  only  once,  and  has  complexity  0(n.3)  [3. 
cap.  4],  □ 

4.  CONCLUDING  REMARKS 

The  algorithm  described  in  the  previous  section  assumes  that  the  target  is  given.  Two  criteria 
seem  to  be  suitable  to  select  a  "good"  target: 

1  1  >  maximum  covering  of  the  dead  cells; 

(  2  i  uniform  distribution  of  the  live  cells  outside  the  target. 
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Criterium  1  is  preferable  if  the  number  of  uncovered  cells  is  o  (  I  D  i ). 

A  variation  of  the  problem  arises  when  the  number  of  live  cells  is  not  a  perfect  square  and  all 
the  live  cells  have  to  be  used  in  the  reconfiguration  process.  In  this  case  the  reconstructed  array  is  a 
square  array  that  misses  some  border  cells,  and  the  corresponding  intermediate  array  can  be  con¬ 
structed  with  a  slight  modification  of  the  above  technique. 

The  network  approach  can  be  extended  in  order  to  take  into  account  also  the  channel  width. 
Let  A  U  "  *' be  the  network  obtained  from  *  by  assigning  capacity  c  j  to  each  node  and  capa¬ 
city  c2  to  each  arc  in  A}k  \  By  using  arguments  similar  to  the  ones  of  section  2.  we  obtain  that  if 
A/(1‘ ,  c  2'  admits  a  flow  of  value  I  D  p)  (V  —  7")  I ,  then  there  is  a  reconstruction  of  the  array  with 
maximum  wire  length  and  channel  width  bounded  by  c  j  and  c2.  respectively. 
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